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ABSTRACT 

Aims. To study the kinematics and dynamics of the extreme outer halo of M87, the central galaxy in the Virgo cluster, and its transition 
to the intracluster light (ICL). 

Methods. We present high resolution FLAMES/VLT spectroscopy of intracluster planetary nebula (PN) candidates, targeting three 
new fields in the Virgo cluster core with surface brightness down to \i B = 28.5. Based on the projected phase space information 
(sky positions and line-of-sight velocities) we separate galaxy and cluster components in the confirmed PN sample. We then use the 
spherical Jeans equation and the total gravitational potential as traced by the X-ray emission to derive the orbital distribution in the 
outer stellar halo of M87. We determine the luminosity-specific PN number for the M87 halo and the ICL from the photometric PN 
catalogs and sampled luminosities, and discuss the origin of the ICL in Virgo based on its measured PN velocities. 
Results. We confirm a further 12 PNs in Virgo, five of which are bound to the halo of M87, and the remainder are true intracluster 
planetary nebulas (ICPNs). The M87 PNs are confined to the extended stellar envelope of M87, within a projected radius of ~ 160 
kpc, while the ICL PNs are scattered across the whole surveyed region between M87 and M86, supporting a truncation of M87's 
luminous outer halo at a 2cr level. The line-of-sight velocity distribution of the M87 PNs at projected radii of 60 kpc and 144 kpc 
shows (i) no evidence for rotation of the halo along the photometric major axis, and (ii) that the velocity dispersion decreases in the 
outer halo, down to <x/ fl « = 78 + 25 km s _1 at 144 kpc. The Jeans model for the M87 halo stars fits the observed line-of-sight velocity 
dispersion profile only if the stellar orbits are strongly radially anisotropic (J3 ^ 0.4 at r — 10 kpc increasing to 0.8 at the outer edge), 
and if additionally the stellar halo is truncated at =i 150 kpc average elliptical radius. The a-parameters for the M87 halo and the ICL 
are in the range of values observed for old (> 10 Gyr) stellar populations. 

Conclusions. Both the spatial segregation of the PNs at the systemic velocity of M87 and the dynamical model support that the stellar 
halo of M87 ends at ~ 150 kpc. We discuss several possible explanations for the origin of this truncation but are unable to discriminate 
between them: tidal truncation following an earlier encounter of M87 with another mass concentration in the Virgo core, possibly 
around M84, early AGN feedback effects, and adiabatic contraction due to the cluster dark matter collapsing onto M87. From the 
spatial and velocity distribution of the ICPNs we infer that M87 and M86 are falling towards each other and that we may be observing 
them just before the first close pass. The new PN data support the view that the core of the Virgo cluster is not yet virialized but is in 
an ongoing state of assembly, and that massive elliptical galaxies are important contributors to the ICL in the Virgo cluster. 

Key words, galaxies: clusters: individual (Virgo) — stellar dynamics — (ISM:) planetary nebulae: general — galaxies: halos — 
galaxies: elliptical and lenticular, cD — galaxies: formation 



1. Introduction 



Over the past few years the diffuse intracluste r light (ICL) has 
been t he focus of many st udies, both in nearby dFeldmeier et all 



2005 ) and in intermedia te redshift clusters 



2004 ; iMihoset all 

dZibetti et alll2005l;TKrick & Bernsteirll2007l) . It has been found 
that the ICL is centrally concentrated and in many cases, includ- 
ing the diffuse outer halos o f galaxies, comprise s ~10% of the 
total starlight in the cluster ( Zibetti et all I2005I) . and up to as 
much as ~ 35% dGonzalez et alll2007l) . 



Theoretical studies of the diffuse cluster light through sim- 
ulations predict that the ICL is unmix ed and therefore should 
exhibit a fair amount of sub-struc t ure ( Napolita no et al 1 12001 
iMurante et all 120041: iRudick et all 120061) . An important contri- 
bution to the diffuse light in clusters may come from the ex- 
tended halos of giant galaxies: numerical simulations predict 
the presence of such halos around isolated galaxies out to sev- 
eral hundred kpc, c onsisting of stars shed by merging sub-units 
(lAbadi et all 120 06). When these galaxies enter the cluster core, 
their halos would be stripped first by the tidal fields and later 
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by the tidal shockin g in the interaction with the cluster's centr al 
core and cD galaxy (R udick et ail 120061: iMurante et alll2007l) . 

Indeed t he de ep image of the Virgo cluster core by 
Mih os et ail (|2005), reaching yuy = 28 mag arcsec~ 2 , shows a 
variety of features such as streamers, arcs and smaller features 
associated with individual galaxies. It also shows faint, very ex- 
tended diffuse halos surrounding the large g alaxies. In par t icular , 
around the giant elliptical galaxy M87, the [Mihos et all d2005l) 
photometry reveals an extended stellar envelope at very low sur- 
face brightness levels, fiy > 26.5 mag arcsec~ 2 , with flattened 
isoph otes (noted previously by Arp & Bertolal,[i97lHWeil et all 
1 1 997b . and out to ~ 37' 161 kpc) along the semi-major axis. 

The Virgo cluster has long been known to be dynamically 
unmixed, with complex sub-structures. This was first realized 
from the spatial and velocit y distribution of Virgo galaxies (e.g. 
iBinggeli et all 1 19871 [l993b . In particular Binggeli et al. (1993) 
found tentative evidence from the asymmetry in the velocity dis- 
tribution of dwarf spheroidal galaxies that even the core of Virgo 
is not virialised, and suggested that the cluster is dynamically 
young, with two sub-clumps M87 and M86 falling in towards 
each other in the centre. 

From photometry in z ~ 0.1 clusters dGonzalez et all 12005; 
iRrick & Bernsteirl | 2007l) and from kinematic studies of the ICL 
in nearby clusters dArnaboldi et all 12004 iGerhard et ail I2005L 
120071) we have learned that the genuine ICL component, defined 
as the light radiated by stars floating freely in the cluster poten- 
tial, and the extended halos of bright (elliptical) galaxies often 
overlap spatially, and cannot easily be distinguished from broad- 
band photometry alone. Kinematic information can complement 
the photometry. For surface brightness hb < 25mag/arcsec 2 , in- 
tegrated light spectroscopy can be used to measure the mean ve- 
locity and veloci ty dispersion in the outer halos of the brightes t 
cluster galaxies dSembach & Tonrvl 119961; iKelson et all 120021) : 
however, reaching the faint surface brightness level of the true 
ICL component with this technique is very difficult. Since plan- 
etary nebulae (PNs) follow light (e.g. lCoccato et all [20091) . the 
spectroscopic study of these tracers, both in the extended halos 
and the ICL, offers a way to identify and measure the kinematics 
of these diffuse stellar components down to very faint surface 
brightness (fig < 28.5 in Virgo), but it is currently limited to 
clusters with distances < 100 Mpc ( Gerhard et alll2005l) . 

For the Virgo cluster, there has been considerable success 
with a two-step approach of identifying PN candidates with 
narrow-band imaging followed by multi-object spectroscopy. 
lArnaboldi et al.l d 19961) observed the outer regions of the giant el- 
liptical M86, measuring velocities for 19 objects. Three of these 
turned out to be true ICPNs, with velocities similar to that of 
the mean velocity of the Virgo cluster. Subsequently, 23 PNs 
were detected in a spectroscopic survey with 2dF on the 4m 
Anglo-Australian Telescope (Freeman et al. 2000; Arnaboldi et 
al. 2002). These results were all based on single line identifica- 
tions, although the second oxygen line was seen with the right 
ratio in the composite spectrum of 23 PNs observed by Freeman 
et al. (2000). The first confirmation based on dete cting the TOUI! 
doublet in a single PN spectrum was made in lArnaboldi et alj 
(2003). Expanding on this early work, we began a campaign to 
systematically survey PN candidates in the Virgo cluster using 
multi-object spectrosco py with the FLA MES/GIRAFFE spec- 
trograph on the VLT dArnaboldi et all 12004 hereafter A04). 
A04 presented the first measurements of the velocity distribu- 
tion of PNs from three survey fields in the Virgo cluster core and 
concluded that in two of these fields the light is dominated by 
the extended halos of the nearby giant elliptical galaxies, while 
the ICL component dominates the diffuse light in only one field, 



where a 'broad' line-of-sight velocity distribution is measured, 
and all PNs are true 'ICPNs' . 

We here present PN velocity measurements from a further 
three pointings in the heart of the cluster core. We emphasize 
that these pointings are targeting faint surface brightness regions 
well outside of individual galaxies, in order to trace the ICPNs 
expected to be moving freely in the cluster potential, and thus to 
investigate the dynamical state of the ICL and of the core of the 
Virgo cluster. The photometric/geometric classification of PNs 
as ICPNs is in fact revised later in the paper according to the dy- 
namical information obtained from the line-of-sight (LOS) ve- 
locities of the confirmed PNs. 

In this paper we give a summary of our observations and 
data reduction in § [2] where we also discuss the sample com- 
pleteness and show the final emission spectra. The distribution 
of measured line-of-sight velocities (LOSVD) and the projected 
phase-space diagram for these PNs are presented in § [3] From 
these data we distinguish between ICPNs and PNs bound to the 
halo of M87. In §[4] we discuss the rotation, velocity dispersion 
and physical extent of the stellar halo of M87, using the veloci- 
ties of the PNs bound to M87 in the combined data sets of this 
paper and A04. We then construct a dynamical model based 
on the gravitational potential obtained from X-ray observations 
and the combined absorption-line and PN velocity dispersion 
data for the galaxy. Possible mechanisms for the truncation of 
M87's stellar halo are discussed in § 14.51 In §[5]we compute the 
luminosity-specific PN number 0-2.5 for both the M87 halo and 
the ICL in Virgo, and in § [6] we discuss the implications of the 
ICPN LOSVD for our understanding of the dynamical status of 
the cluster core and the origin of the ICL in Virgo. Summary and 
conclusions of the paper are given in § [7] 

In what follows, we adopt a distance of 15 Mpc for M87, 
equivalent to a distance modulus of 30.88; therefore 1" =73 pc. 

2. Observations 

The observations were taken in service mode (22 hrs, 076. B- 
0086 PI: M. Arnaboldi) over the nights 25-28th March 2006 us- 
ing the FLAMES spectrograph on UT2/VLT in MEDUSA mode 
which allows spectra to be taken through up to 132 fibers simul- 
taneously The data were taken in clear conditions with seeing 
< 0.9". We used the high resolution grism HR504.8 centred at 
504. 8nm and with wavelength coverage 250A and spectral res- 
olution 20,000. With this setup, the instrumental broadening of 
the arc lines is FWHM = 0.29 A or 17 kms 1 , and th e error on the 
wave length measurements is 0.0025 A or 150 ms _1 dRover et all 
120021) . 

Figure 1 shows the location of the selected fields targeted 
with FLAMES, including the three previous fields FCJ, CORE 
and SUB presented in A04, and the three new fields F4, F7_l, 
F7_2. The photometry used fo r the selection of PN candidates 
is from Feldmeier et al. (2003); his fields F4 and F7 contain the 
FLAMES fields F4 and (F7.1, F7_2), respectively. 

2.1. Data reduction and sample completeness 

The data were reduced using the GIRAFFE pipelin^] including 
bias subtraction, determining fiber location on the CCD, wave- 
length calibration, geometric distortion corrections and extrac- 

1 See http://www.eso.org/sci/facilities/paranal/instruments/flames/ 
overview.htmT. 

2 The GIRAFFE pipeline is available at 
http://girbldrs.sourceforge.net 
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Fig. 1. Deep image of the Virgo cluster core showing the diffuse light distribution (Mih os et al.U 2005t). with our target fields super- 
posed. Target fields of the previous spectroscopy (A04) are shown as red circles and our new target fields as well. The blue ellipse 
shows the boundary used in the dynamical modeling in § 14.41 



tion of the one-dimensional spectra. The co-addition was carried 
out separately as a final step on the one-dimensional spectra as 
the fibers are allocated in a slightly different order for MEDUSA 
plates 1 and 2 and the pipeline does not account for this. 

Table 1. Observed Fields and Spectroscopic confirmations 



Field 


F7.1 


F7_2 


F4 


a(J2000.0) 


12 28 53.70 


12 28 46.40 


12 27 43.35 


<S(J2000.0) 


+ 12 44 32.5 


+ 13 00 20.5 


+ 12 33 57.5 


ly tot 


22 


28 


25 


N* 


13 


12 


8 


N c 

1 1 em 


9 


5 


5 


N'' 


6 


4 e 


3 



" Total number of targets with allocated fibers. 

Number of targets with allocated fibers and m J0 07 brighter or equal to 
the completeness magnitude limit of the photometric survey. 
c Number of spectra with detected emission line. 
d Number of spectra with both [O m]A 4959, 5007 A detected. 
e One of these is in common with F7_l, i.e., we have only 12 confirmed 
PNs in total, one of which was observed twice. 




2b.b 26.G 27.G 27.S 



Fig. 2. Histogram showing the msoo7 magnitudes of all our ob- 
served targets (solid black line) over-plotted with those where 
emission lines were detected (red dashed line). The blue dotted 
lines show the photometric completeness limits for target fields 
F4 (26.6) and F7 (26.8). 



Table Q] shows the number of spectroscopically confirmed 
emission-line objects and planetary nebulas, with respect to the 
number of candidates targeted (N, ol ), and the number of candi- 
dates targeted above the photometric completeness limit for each 
field (Ncmp). 

A histogram showing the number of candidates versus num- 
ber of confirmed emission-line objects by magnitude is shown 
in Figure [2] The photometric completeness limits (~ 90%) in 



the two photometric fields F4 (m5007=26.6) and F7 (mgoo7=26.8) 
are shown as blue dotted lines. These photometric complete- 
ness magnitude limits are defined as to where the signal-to-noise 
over the entire photometric measurement is nine per pixel or 
greater, correspo nding to a photometric error of approximately 
0.12 magnitudes dFeldmeier et al.Ll2003l) . The confirmation rate 
for emission-line objects above the completeness limit is then 
40-70% depending on the field. This is comparable with the re- 
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Fig. 3. Spectra for the confirmed PNs, ranked by magnitude 1115007 and smoothed by a factor 7 to 0.035nm per pixel. 1115007 and 
the LOS velocities are labelled in the top left corner of each spectrum. The expected location of [Om]/l4959 is shown by a red 
dash, in the cases where not immediately obvious. The [Oiii]5007 line for the PNs with m.5007=27 and v/ os = -155km s has a 
signal-to-noise ratio of 3 and this is our weakest believable detection. A typical error in the velocities is 3 km/s; see Section |Z21 



suits from A04 (30-80% varying by field), and is a reasonable 
recovery rate given the following effects. 

Firstly, if the astrometry is not very precise or if some rota- 
tion error is introduced in positioning the plate, then part or all of 
the flux from some objects may miss the fibers. This is more se- 
rious for faint objects as they will then not be detected above the 
noise. Indeed, in Figure[3]the total flux in the [Om]5007 A line 
is not clearly correlated with the magnitude 1115007 of the source, 
indicating that fiber positioning might be problematic. We mea- 
sure the relative fluxes for the same [Om]5007 A detection in 
different frames and find that it can in fact vary by up to a factor 
of two. 

The likelihood of having false candidates above the com- 
pleteness limit is very low. Each candidate ICPN was hand- 
inspected, and the code that finds the objects has been ex- 
tensively tested on closer galaxies. In some cases (e.g., M51, 
iDurrell et alJ d2003)). there has been close to 100% recovery, us- 
ing the same techniques as used to select the candidate ICPNs 
here. 

However, below the photometric completeness limit the un- 
certainties are clearly much higher. Although many of the fainter 
objects are still likely to be ICPNs and were hence targeted 
spectroscopically, the probability for 'spillover' (Aguerri et al. 



2005) increases substantially. Due to the photometric errors in 
the [OIII] fluxes some objects will be measured with a brighter 
flux than their real flux. If in addition their broad-band fluxes fall 
below the limiting magnitude of the off-band image they will be 
selected as ICPNs when they are in fact very faint continuum 
stars, due to the fact that they will have erroneously large nega- 
tive colours. 



2.2. Spectra of PNs and background emission line galaxies 

Figure [3] shows all of the spectra for the confirmed PNs, 
ranked by their photometric magnitude 015007- F° r most of the 
PNs brighter than 1115007=27 we also detect the second line 
[Om]/l4959. The expected location of [Om]/i4959 is shown by 
a red dash, where not immediately obvious by eye. The target 
fields F7_l and F7_2 overlap and have one source in common. 
The independently measured velocities for this source agree to 
within 3km s _1 . 

Our weakest believable detection has a total signal-to-noise 
ratio S/N=3. As an additional check we create the average com- 
bined spectrum for the 12 identified PNs (Figure |4| and mea- 
sured the equivalent width ratio of the two [OIII] lines. The ratio 
is 3 as expected if all identifications are real. 
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Fig. 4. Combined spectrum of all 12 identified PNs, Doppler cor- 
rected to the rest-frame. 
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Fig. 5. Examples of the other emission-line objects present in 
the sample, [Oil], Ly— a, and an unidentified broad emission line 
which we speculate might be Hell at 1640 A in a high-z LBG, or 
alternatively CIV 1550 or [CIII] 1909 from an AGN. The spectra 
have been smoothed to 0.035 nm per pixel. For comparison to 
Figure [3] the 015007 magnitudes are shown in the top right corner. 



In Figure [5] examples of the other emission-line objects 
present in the sample are shown: an [0 11] doublet, an asymmetric 
Ly-a line, and an unidentified broad emission line which might 
be an AGN (for example CIV1550 or CIII] 1909). Alternatively, 
there is a possibility that the line is Hell at 1640 A in a LBG at 
higher redshift. Shapley et al. (2003) discuss that this is some- 
times seen as nebular emission and also as much broader emis- 
sion (~ 1500km s -1 ) possibly from stellar winds. The FWHM 
of the lines we see is about ~ 200km s (observed frame). The 
contamination rate of these other emission-line objects is 7/19, 
or 37%. 



500 1000 1500 

Velocity (km/s) 



Fig. 6. Radial velocity histograms. The bottom panel shows the 
velocity distribution of all identified PNs in the 3 new fields. 
The peak at ~ 1300km s _I corresponds to PNs bound to the 
halo of M87. They have a mean velocity of 1297 + 35km s 
and rms dispersion of 78 + 25km s . The middle panel shows 
the distribution of PN velocities in the previously surveyed FCJ 
field (A04), and the top panel shows the distribution of dwarf 
spheroidal LOS velocities in the same region of the Virgo clus- 
ter core, for comparison. The systemic velocities of M87 and 
M86 are shown with the dotted lines. 



The PN emission lines are all resolved, and thus we have 
been able to measure the expansion velocities of these PNs and 
to derive information on the masses of t he progenitor stars. Thi s 
work is presented in a companion paper (Arnab oTdi et all [2008 ). 
Here we concentrate on the kinematics, yielding information on 
the halo of M87 and the assembly history of the Virgo cluster. 



3. LOSVD and projected phase-space 

Figure[6]shows the distribution of velocities of the newly identi- 
fied PNs in the Virgo cluster core. The velocities have not been 
adjusted for a heliocentric correction, as the observations were 
almost all taken close to the equinox and the correction is within 
only +3km s _1 . 

For the subsequent analysis, we combine these velocities 
with the A04 sample in the FCJ and Core fields. Figure [7] shows 
the l ocation of these PN s on the deep image of the cluster core 
from Mih os et all (120051) . and Figure [8] shows their distribution 
in the projected phase-space plane defined by projected distance 
from M87 center and line-of-sight velocity 

In the phase-space diagram Figure|8] we can identify two re- 
gions with very different characteristics: For projected distances 
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Fig. 7. Deep image of the Virgo cluster core showing the distribution of the intracluster light (Mihos et al. (2005)). The spatial 
distribution of our spectroscopically confirmed PNs are overlaid. The A04 targets are shown in green. Our new targets are shown in 
red if redshifted with respect to Earth and blue if blueshifted. Objects with velocities higher than the mean velocity of Virgo (1 100 
km s _1 ) are shown as crosses and those with lower velocities shown as circles. Dwarf spheroidals are marked as magenta dots. The 
velocities (in km s" 1 ) are labelled for all objects shown. The nominal 'edge' of the M87 halo at ~~ 160 kpc is indicated with a 
white circle. The pink circle has a 1.5 degree diameter and is centered on the projected midpoint of M87 and M86. North is up and 
East is to the left. 



R < 2400" most of the PNs are strongly clustered around the 
systemic velocity of M87, v jyj = 1307 km s . By contrast, for 
R > 2400", the PN velocities spread widely over a velocity range 
more typical for the Virgo cluster. From the latter, intracluster 
region we see a string of low PN velocities (800-400 km s -1 ), 
extending inwards to the upper FCJ field (see Fig.|7]i. 

In the FCJ field at projected distance R < 1300" there are 
two of these intracluster outliers at ~ 400 km s . The remain- 
ing PNs are distributed symmetrically around v SVJ and have mean 
velocity 1276 + 71km s and velocity dispersion <x = 247 
km s _1 (A04); their velocity distribution is shown in the middle 
panel of Fig. [6] 

In the combined new F7/F4 fields at R ~ 2000" we find 
five PNs tightly clustered around v JVJ = 1307 km s _1 ; these 
have mean velocity 1297 ± 35km s and an rms dispersion of 
78 + 25km s . At comparable radii there are two additional PNs 
with velocities of 753 and 634 km s _1 ; compared to the previous 
five, these are Icr and 8cr outliers. It is unlikely that one or two of 
these outliers are part of the same (very asymmetric) distribution 
as the five PNs clustered around v sys . By contrast, they fit natu- 
rally into the stream leading from the FCJ outliers all the way 
into the ICL. We therefore identify as PNs bound to the M87 
halo only those PNs which are clustered around the systemic ve- 
locity of M87. These are confined to radii R < 2400". 



The M87-bound PNs in the FCJ and combined F7/F4 fields 
are located at mean projected radii of 60 and 144 kpc, respec- 
tively. They correspond to the narrow peaks in the line-of-sight 
velocity distributions (LOSVD) in the lower and middle his- 
tograms in Fig. [6] 

Outside R = 2400" in Fig. [8] we find PNs at larger relative 
velocities to M87, with an approximately uniform distribution 
in the range -300 to 2600 km s _1 . Those in the radial range 
2400" < R < 3600" (the F7/F4 field) are confined to negative 
velocities with respect to M87. These are probably encroaching 
stars from M86 and other Virgo component^ By contrast, the 
PNs further than 3600" from M87 (in the Core field) show a 
broad distribution of velocities, more characteristic of the cluster 
as a whole (see A04). 

In the middle and bottom panels of Fig. [6j the ICL PNs 
show up as approximately flat velocity distributions in their ve- 
locity range, besides the peak of velocities from PNs bound to 
M87. A flat distribution of velocities in addition to the peak 
near M87's systemic velocity is also seen in the LOSVD of the 
dwar f spheroidal galaxies i n the same region of the Virgo cluster 
core <BinggelietaflLl993|> which is shown in the top panel of 
Figure[6] However, for the dwarf galaxies the flat velocity distri- 

3 They cannot be in the Local Group: th e faintest PNs in the SMC 
have ni5oo7=23 I Jacob v & De Marcel [2002) so are still much brighter 
than the brightest Virgo PNs at m J00 7=26.3. 
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Fig. 8. Distribution of line-of-sight velocity versus projected 
distance from the center of M87 for all spectroscopically con- 
firmed PNs in the new fields as well as the FCJ and Core fields 
of A04. 

bution extends to significantly more redshifted velocities, indi- 
cating that the dwarfs and ICL PNs kinematics can only partially 
be related. 



4. The M87 Halo 

We have seen from the phase-space diagram in Fig. [8] that the 
PNs in the FCJ and F7/F4 fields divide into two components, one 
associated with the halo of M87, and the other with the unbound 
Virgo ICL. All PNs found around the M87 systemic velocity are 
within R = 161 kpc projected radius from the galaxy's center. 
In the following subsections we combine our velocity measure- 
ments with the kinematic data in the literature, and discuss the 
rotation, velocity dispersion profile, dynamics and truncation of 
the outer M87 halo. Finally we consider possible origins of the 
truncation. 



4.1. Rotation of outer halo? 

First we ask whether there is any evidence in our data for rotation 
in the out er halo of M87. For the globular cluster (GC) system 
of M87, ICohen & Rvzhovl dl997l) inferred a rotation of about 
100 km s _1 for R < 35 kpc, approximately about the minor axis 
of the galaxy intensity isophotes, using spectra of low resolution 
with errors for the GC velocities of order 100 km s _1 . Cohen 
(2000) added new data for GCs in the halo at 24 < R < 43 kpc 
with smaller errors (typically ~ 50k m s _1 ), and inferred a rota- 
tion of 300km s~ 1 . ICote et al J d200 ll) carried out an independent 
analysis using a ne w spectroscopic and photo metric database 
dHanes et all |2001[) partly based on that of ICohen & Ryzhov 
j 19971) an d (Cohenl d2000l) . and similarly find ~ 160km s" 1 . 
ICote et all d2001l) found that the metal rich GCs rotate every- 
where about the photometric minor axis of the galaxy, while 
the metal poor GCs have a more complex behaviour: they ro- 



tate about the photometric minor axis of the galaxy between 
15 < R < 40 kpc, and about the major axis at R < 15 kpc. 

If the PN population in the outer halo of M87 also rotated 
about the galaxy's photometric minor axis, similarly to the M87 
GC system at 15 < R < 40 kpc, we should see a signature along 
the major axis of the galaxy, that is, the mean velocities of the 
PN LOS VD peaks a ssociated with M87 in our two pointings FCJ 
(A04; this is F3 in lFeldmeier et aUl2003l) and F7/F4 (Figure 0) 
should be offset from the systemic vel ocity of the galaxy (1307 



km s '). From the extrapolated fit of ICote et alj d200ll) to the 



GC radial velocities we would expect a negative constant offset 
of about 160 km s~'at both field positions, i.e., a mean velocity 
of ~ 1150 km s _1 . 

For the M87 sample of PNs identified in the phase-space dia- 
gram in Fig. [8] we find a mean velocity of v ra d = 1297+35km s 
in the new field (F7/F4) at mean projected radius R = 144 kpc. 
In the previously studied field (FCJ; A04) the mean velocity is 
1276 + 71km s" 1 at mean R — 60 kpc. Thus we see no evidence 
for rotation of the outer stellar halo around the galaxy's minor 
axis in either the PN sample at R = 60 kpc or at R = 144 kpc. 
The rotation seen in the GCs may thus suggest that they do not 
trace the main stellar population of M87, or that they are con- 
taminated with IC GCs with a LOSVD similar to the encroach- 
ing st ars of M86 as seen in Fig.[8](see also Fig. 1 of ICote et all 
2001). We have not surveyed fields along the minor axis, so we 
cannot check with PNs whether there is rotation about the major 
axis. 



4.2. Velocity dispersion profile 

With the new data, we can also now plot the velocity dispersion 
profile of M87 all the way out to R = 144 kpc from the centre 
of the galaxy. Figure [9] shows this as a function of projected ra- 
dius from M87 centre. In the inner regions (R < 25") we use the 
G-band absorption line measurements from the integrated stellar 
light of van der Marel (1994). In the region 25" < R < 80" we 
use stellar velocity dispersions from Sembach & Tonry (1996). 
As these authors discuss, their data is systematically offset from 
most other datasets by 7 - 10%, d ue to using a larger slitwidth. 
Romanowsky & Kochanek (2001) calculate that this amounts to 
an additional instrumental dispersion of 183km s _I and so we 
adjust the velocity dispersion by this amount (subtracting in 
quadrature) to bring it in line with the van der Marel dataset. We 
take the average of the velocity dispersions at each positive and 
negative R, assuming symmetry with respect to the galaxy's cen- 
ter. We exclude the Sembach & Tonry data in radial bins beyond 
80" as there is a discrepancy between the velocity dispersions at 
the corresponding opposite positions in radius along the axis and 
furthermore the values in those bins have large error bars. This 
may be due to low S/N in the outer part of the galaxy where the 
surface brightness is low and/or real anisotropics in the velocity 
dispersions. Either way we judge it better to exclude these data 
points as they are less trustworthy. 

We also show in Figure|9]the data of lCote et all d200ll) for the 
metal-rich GCs out to 380". We exclude the outer bins (380"- 
635") where the error bars are close to 100% and therefore do 
not constrain the shape of the velocity dispersion profile in any 
way. The metal-poor GC system is more spatially ext ended than 
its me tal-rich counterpart (see Figure [10] below and ICote et all 
120011) and may be composed of accreted and/or infalling remains 
of 'failed' or disrupted dwarfs: their velocity dispersion data are 
also shown in Figure [9] but do not trace the veloc ity disper- 
sions o f the M87 stars. We also note from Figure 1 of Cot e et all 
(2001) that the GC sample is likely to contain intrac luster GCs 
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Fig. 9. Velocity dispersion profile of M87, including stellar ve- 
locity dispersions from absorption-line spectra and discrete LOS 
velocity measurements from glob ular cluster and PN data. The 
squares are data poin ts fromlvan der Marell (Il994l) . the green di- 
amonds are based on lSembach & Tonrvl d!996l) . and the red and 
blue stars are velocit y dispersions for th e metal-rich and metal- 
poor GC samples of ICote et alj d200ll) . respectively. The ma- 
genta circle is the PN velocity dispersion determined in A04 and 
the magenta triangles are the PN dispersion values from this pa- 
per. These last two points are approximately along the major axis 
of the outer isophotes, which have ellipticity e - 0.43. 



just as our PN sample contains ICPNs, requiring a careful anal- 
ysis of the GC LOS VDs. 

Finally, the two outermost velocity dispersion points are 
from planetary nebulas presented in A04 and this paper. We note 
that when we bin the A04 data to be consistent with the bin- 
ning of the velocity distribution in this work (100km s bins) 
the peak around M87 is resolved into a somewhat narrower peak 
of 9 objects, with two lower velocity and one higher velocity 
outliers (Figure [9). The mean and rms velocity of this peak of 
v = 1264 km s _I and cr = 247 km s^'from A04 then change to 
1292 + 46km s and 139 + 23km s respectively. It is possible 
that the larger value of ctpn obtained by A04 could be due to 
the inclusion of some ICPNs from the component with uniform 
LOSVD seen in Figure|6] We carried out x 2 testflbut could not 
distinguish between both interpretations. The dynamical model 
discussed below favors the high value of ctpn at R — 60 kpc. 

In any case, the PN velocity measurements show that the 
halo of M87 becomes colder at larger radii: the velocity disper- 
sion decreases to 78km s _I at R ~ 140 kpc. 



4 We carried out a x 2 test f° r me (FCJ; A04) sample and i) a broad 
Gaussian (v = 1264 km s and cr = 247 km s" 1 ), ii) a uniform distri- 
bution plus a narrow Gaussian (v = 1264 km s and cr = 139 km s _I ). 
Because of the limited statistics of the PN sample in this field, the re- 
sults depend on the velocity range chosen for the test. In a 700 - 1650 
km s - ' range, the two distributions fit the data equally well with 80% 
probability, while the broad Gaussian is ruled out in a 350- 1650 km s" 1 
range. 
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Fig. 10. V-band surfa ce brightness profile for M87 from 
Kormend v et al.l (J2008) shown with black circles along aver- 
age ellipse radii R avg of the isophotes. Red circles show the ex- 
trapolated Sersic fit. The full line shows the reprojected three- 
dimensional luminosity model obtained from the extrapolated 
data, and the dashed line shows the reprojection of the luminos- 
ity model when truncated at average ellipse radius R trU nc = 149 
kpc. For comparison, the number dens ity profiles of the red and 
blue globular cluster populations from lCote et alj d200ll) are also 
shown with arbitrary scaling as red and blue stars, respectively. 



4.3. Truncation of the M87 stellar halo 

In the FCJ field, there are M87 halo PNs detected throughout, 
but in the F7/F4 fields, the PNs around the systemic velocity 
of M87 (=1307 km s _1 ) appear to be found only within a pro- 
jected radius of R = 161 kpc (see the spatial distribution of the 
spectroscopically confirmed PNs in Figure|7]). At projected radii 
R > 161 kpc from the center of M87, we find only the encroach- 
ing stars of M86 and other Virgo components. We now investi- 
gate whether this spatial and velocity segregation is significant 
and indicates that the M87 stellar halo is truncated. 

iKormendv et alj d200 8) present a composite V-band surface- 
brightness profile for M87 out to 135 kpc along the semi-major- 
axis. This is shown with black circles in Figure [TUJ with the 
semi-major axis of each isophote replaced by the average el- 
lipse radius, R avg = (ab) 1 ^ 2 of the isophote. We will use the 
latter in the constructio n of the spherical Jeans models in § 14.41 
IKormendv et alj d20 08) fit a Sersic profile to the semi-major axis 
profile excluding the central core and the last two data points 
(which may have a significant ICL contamination), and obtain 
the following best-fit parameters: /i c = 25.71, R e = 704" = 51.2 
kpc, n = 11.885. We use this Sersic fit to the surface bright- 
ness profile (see Fig. |T0b to compute the luminosity of the M87 
halo at radii outside the available photometry. We note that the 
description of M87 as a classical E0 or El galaxy is based on 
short exposure optical images, while in deep images its isophotes 
show marked eccentricity. For the extrapolation we assume an 
ellipticity e = 0.43 and position angle (PA= -25°), based 
on the outer parts of the ellipticity and position angle profiles 
in IKormendv et al.1 d2008l) . We then evaluate the M87 halo lu- 
minosities in the regions of overlap between the photometric 
and spectroscopic fields, in which, respectively, the photometric 
identification and spectroscopic follow-up of the PNs was car- 
ried out. These are shown by the colored regions in Figures QT| 
For ellipticity e = 0.43, the PNs belonging to M87 appear to be 
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Fig. 11. Left: The red region is the intersection between the pho- 
tometric field FCJ and the FLAMES (FCJ) pointing. Right: The 
red region is the part of the intersection between the photomet- 
ric field F7 from Feldmeier et al. (2003) and the regions jointly 
covered by the F7_l, F7_2 and F4 FLAMES circular pointings, 
which is within the isophote with an average ellipse radius, 
Rtmnc - 149 kpc. This covers the region containing PNs with 
velocities bound to M87. The green region is the part of the in- 
tersection between the F7 region and the FLAMES pointings as 
before, but outside R m mc, where no PNs associated with the M87 
velocity peak are found. 

found within an average ellipse radius R avg = 149 kpc. We use 
the isophote corresponding to this radius to demarcate the region 
containing PNs belonging to M87 (in red in Fig.fTTT) from the re- 
gion containing no PNs belonging to M87 (in green in Fig.fTTI). 
Table [2] gives the areas and the V-band luminosities of the vari- 
ous regions obtained. On the basis of the Sersic fit, the ratio of 
the M87 luminosity in the F7-green and in the F7-red area is 
0.92. 

From the number of spectroscopically confirmed M87 PNs 
in the F7-red area (i.e., 5), and the ratio of the M87 halo lumi- 
nosities from the Sersic fit (i.e., 0.92), we can then predict the 
expected number of PNs at the M87 systemic velocity in the F7- 
green area, to be 5 x 0.92 = 4.6 PNs. The observational result of 
zero M87 PNs detected in the F7-green area thus implies a trun- 
cation of the M87 stellar halo at R avg > 149 kpc, at a ^ 2cr level. 
This radius tells us the location of the outermost PN in terms of 
the average ellipse radii and we now refer to it as the truncation 
radius R trunc . 

Table 2. M87 halo PNs and sampled luminosities. 



Field 


^spectr 


Area 


L 




kpc 2 


Lv,s 


FCJ red 


12 


4157 


8.7 x 10 9 


F7 red 


5 


4349 


1.3 x 10" 


F7 green 





8909 


1.2 x 10" 



This is a surprising result: the numerical simulations of 
lAbadi et alj (12000) find that the luminous halos around isolated 
galaxies should extend to the virial radius, i.e., to several hun- 
dred kpc and well beyond their traditional luminous radius. It 
is therefore not obvious why there should be a truncation of 
the M87 stellar halo, see Section 14.51 One question is whether 
the truncation we see occurs only at the targeted P. A. within an 
opening angle 6P.A., say, rather than at all azimuths, and whether 
when azimuthally averaged, the halo light distribution would ex- 



tend to larger radii. This could be the case if we had reached 
the radii where the stellar halo of M87 has a signific ant amount 
of sub structure, similar to the outer Milky Way halo (iBell et all 
2008). For example, we might explain the spatial segregation of 
the M87 halo PNs in terms of a cold stellar shell at our field po- 
sition, followed by a steeper surface brightness profile at those 
P. A. whereas at other P.A. the profile would be shallower. This 
would also explain the small <tpn that we measure in our fields 
at R — 144 kpc, as the stars populating shells are near to the 
apocenters of their elongated orbits. 

To assess this we must reconsider the photometric structure 
of M87. Within a semi-major axis of ~ 160 kpc and fiy < 27.5, 
the surface brightness distribution around M87 is well approxi- 
mated by an extended envelo pe with c/a 0.57 dMihos et all 
2005; Ko rmendv et all 120081) . This ellipsoidal component in- 
cludes the diffuse fan of stellar material, which extends along 
the projected southeast semi-major ax is out to ~ 100 kpc 
dArp & Bertolal 1 1 97 it IWeil et all Il997h . but is otherwise fairly 
smooth. At larger radii and fainter surface brightnesses, the light 
distribution is a superposition of the outer halo of M87 and the 
ICL and is brightest in the range of P.A. towards M86 where our 
target fields are. At these radii it does show irregular features 
and some radial streamers, and our fields are large enough to in- 
clude several of these. In fact, some of our outermost M 8 7 PNs 
may be close to an arc-like feature in the iMihos et~aT1 (2005) 
image beyond which litt l e ligh t is seen. However, both the ear- 
lier results of IWeil e t al. ( 1997) who reported the apparent lack 
of sharp-edged fi ne structures around M87, and our independent 
inspection of the IMihos et~a El (120051) image near M87, provide 
no evidence for a large number of "shell-like" features at vari- 
ous azimuths and radii around M87. This is true both inside and 
outside our truncation isophote, and in particular around R ~ 60 
kpc, where the PN data already indicate a falling velocity disper- 
sion profile (see Fig. [9j. 

Further investigation of the extended luminosity distribution 
around M 87 would require q uantitative photometry of the deep 
image of IMihos etafl d2005t) . and a larg e-area and wide-angle 
PN velocity survey to separate the outer halo of the galaxy from 
the ICL with better statistics. 

In what follows, we follow an independent approach and 
test the hypothesis of a truncated stellar halo in M87 dynami- 
cally. We will verify whether we can make a dynamical model 
for M87 in which the stellar velocity dispersion reaches low val- 
ues everywhere around M87, and the total gravitational potential 
is traced by the X-ray emission of the hot gas. 

4.4. The mass distribution and anisotropy in the M87 outer 
halo 

The smooth photometric and X-ray emission profiles indicate 
that the outer halo of M87 is in approximate dynamical equilib- 
rium. With the extended velocity dispersion profile we are now 
able to create dynamical models of M87 to infer the orbital struc- 
ture in the outermost halo. In a spherical system, the intrinsic ve- 
locity dispersions of a population of stars with density j moving 
in a potential <t> are related by the second-order Jeans equation 

j r (j(r)cr 2 r (r)) + -fjcr 2 (r) + j(r)— = (1) 

where the anisotropy parameter fi(r) = 1 - (0"e(r)/cr r (r)) 2 quan- 
tifies the orbital structure of the system. Its value spans a range 
between -oo signifying purely tangential orbits and 1 signifying 
purely radial orbits. These intrinsic quantities are now in terms 
of the three-dimensional radius, r. 
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> 400 



Finally, the radial dependence we adopt for the anisotropy 
profile is given by 



r + r\ 

i 

r + r 2 



(3) 



Nulsen & Boehringer, 1995 (ROSAT) 
Scaled to XMM-Newton 



where r\ < r%, f3 m is the maximum anisotropy reached at large 
r, fi m r\ I r2 is the minimum anisotropy and r-i represents the scale 
of the transition from the minimum anisotropy to the maximum 
anisotropy. This solves the Jeans equation ([TJ as a first-order dif- 
ferential equation with the integrating factor 
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I.F. = j(r)r 2fs >" x (r + r 2 ) 2/! 
where x — r\ /r 2 so that 
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Fig. 12. Circular velocity profiles for M87 from X-ray data. cr r (r) = — J I F ydr' 



(4) 



(5) 



Therefore if we know the potential i>, density of stars j and 
assume a radial dependence for the anisotropy, we can solve for 
the intrinsic radial (07) and tangential {ere) velocity dispersion 
profiles. These can then be projected and compared with the ob- 
served projected velocity dispersion profile to fix the parameters 
of the assumed anisotropy profile. Under the spherical assump- 
tion, the radii of the observables will be the average ellipse radii 
R avg introduced in § 14.31 

One method of deriving the potential of a galaxy is to use 
electron density and temperature profiles obtained from X-ray 
measurements of the thermal bremsstrahlung emission from the 
hot interstellar medium surrounding massive elliptical galaxies. 
If the gas is relatively undisturbed then we can assume that the 
gas is in hydrostatic equilibrium and thus derive the potential. 

iNulsen & Bohringerl (1 19951) use ROSAT data and a maxi- 
mum likelihood method to deduce the most likely mass profile 
in the Virgo cluster core, extending from the centre of M87 out to 
300 kpc. They parametrize this profile with a model composed of 
two (approximate) isothermal mass distributions, one attributed 
to the mass of M87 with a mass per unit length /j = 3.6 x 1O 1O M 
kpc 1 and the other to the dark matter of the cluster with a mass 
per unit length Mo/a - 12.4 x 10 I0 M o kpc -1 and a core radius 
a — 42 kpc. This parametrization is given as 



M(r) = fir + M [(r/a) - arctan(r/a)]. 



(2) 



This mass profile is related to the potential through rd<f>(r)/dr = 
v c (r) 2 /r = GM(r)/r 2 , where v c (r) is the circular velocity at r. It 
is more easily calculated than the potential but also independent 
of distance. The circular velocity profile is shown in Figure Q~2] 
and increases from a minimum circular velocity of 393 kms~' to 
eventually a maximum velocity of 830 kms -1 . For comparison 
the same profile is shown with Mo adjusted so that the maxi- 
mum velocity reached is 700 kms -1 , the value suggested by the 
extrapolation of newer observations taken with XMM-Newton 
dMatsushita et aUl2002l) . 

The density of the stars was obtained through the deprojec- 
tion of the surface-brightness profile of M87. To obtain the in- 
trinsic luminosity density, we adopt the iKormendv et aT] £2008) 
profile (also missing the last two data points), and using their 
Sersic fit, extrapolate to very large radii, see § 14.31 then we em- 
ploy the standard deprojection formula. In Figure [TO] the Sersic 
extrapolation is shown by the red circles and the reprojection of 
the intrinsic luminosity density is shown by the solid black line 
which follows the circles very well. 



Then the projected velocity dispersion, crp, is obtained by pro- 
jecting the intrinsic velocity dispersions along the line of sight. 
These integrals were determined numerically. 

To obtain a dynamical model for the Sersic light distribution 
in the potential implied by Fig. [12] we fixed the minimum and 
maximum anisotropy using constant anisotropy models and then 
employed a x 2 minimization approach to deduce the best-fit r 2 
for the solution in eq. [5] This minimization takes into account 
the long-slit data, the ctpn point at 60 kpc from A04, and the 
new <TpN point at 144 kpc, but not the globular cluster data. The 
best-fit model is shown by the solid black line in Figure [T3] It 
fits the data very well within 6 kpc but it is unable to reproduce 
the low PN velocity dispersions in the outer parts, at R avg - 52 
kpc and R avg = 131 kpc, where it would predict LOSVDs with 
dispersions of 350-400 km s _1 . 

On the basis of the results in § 14.31 we assume now that the 
galaxy's intrinsic luminosity density is truncated at r = R lnmc 
1 49 kpc (i.e., in a spherical system the intrinsic truncation radius 
is the same as the projected truncation radius). The reprojection 
of this truncated intrinsic luminosity density is also shown in 
Figure [K)]by the blue dashed line. We construct a Jeans model 
for the truncated luminosity distribution in the same way as 
above. Finally, to check the influence of the assumed potential 
on the results, we have evaluated one further model, also assum- 
ing a truncation in the luminosity density but using the circular 
velocity that was adjusted to the ana lysis of XMM-Newton ob- 
servations in [Matsushita et all (120021) . The corresponding Jeans 
models are shown in Figure [T3] with the blue and green dashed 
lines respectively, with the second model dipping slightly lower 
in the outer parts, reflecting the lower potential in this region. 

Both truncated models behave as the untruncated Sersic 
model in the center, but fall much more steeply at radii R avg > 30 
kpc, thus being able to reproduce the outermost cr PN data points 
at R avg - 52 and R avg = 131 kpc (corresponding to projected 
radii R = 60 and R = 144 kpc). At R avg - 52 kpc, the truncated 
models predict a velocity dispersion which favours the higher 
value of o- PN , i.e., 247 km s _1 . Figure [l4l shows that the best- 
fit models imply a mildy radially anisotropic orbital distribution 
(/? * 0.15) in the centre that becomes highly radially anisotropic 
in the outer halo {J3 » 0.8). 

We conclude that, under the assumption of spherical sym- 
metry, the Jeans models can only reproduce the low PN velocity 
dispersion measurements in the outer halo of M87 at R avg > 30 
kpc only with a truncation of the intrinsic luminous density. 

In principle, this dynamical argument could be circumvented 
if at the position of our outer fields the stellar halo of M87 was 
strongly flattened along the line-of-sight. In this case o~pn could 
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Fig. 13. Velocity dispersion profiles derived for Jeans models 
with spherical symmetry and surface brightness profiles as in 
Figure [TUJ see text. The velocity dispersion data points are 
shown at their average ellipse radii R avg , computed with the el- 
lipticity profile. For the globular cluster velocity dispersions this 
is not possible; these are not used in constructing the Jeans mod- 
els. Symbols are as in Fig. [9] 
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Fig. 14. Anisotropy profiles for the best-fit models: they imply 
a mildy radially anisotropic orbital distribution (J3 » 0.15) in the 
centre that becomes highly radially anisotropic in the outer halo 
08*0.8). 



be low at these radii independent of a truncation. However, at 
~ 150 kpc radius such a flattening is likely to be local and would 
have arisen from the geometry of accretion, rather than signify- 
ing an angular momentum supported global structure collapsed 
from even A~ l ^20 times further out. Thus the well-mixed, 
three-dimensional stellar halo of M87 would then have ended 
at even smaller radii. Also note that in this case we could still 
not explain the lack of PNs at the M87 systemic velocities for 
radii greater than R, nmc (see Fig.|8]and the discussion in § 14.31 ) 



4.5. On the possible origin of a truncated stellar halo in M87 

Summarizing the last two sections, there are two independent 
and mutually consistent pieces of evidence that the stellar halo 
of M87 ends at R tnmc ^ 150 kpc: the lack of PNs around the 
systemic velocity of M87 beyond this radius, and the very low 
velocity dispersion in the outer halo. 

In well-mixed, dense galaxy clusters it is expec ted that 
galax i es are tidally t r uncate d by the cluster's tidal field dMerrittl 
1984; iGhigna et all 119981) . The tidal effects are strongest 
near the cluster core radius and the approximation r t id a i - 
TperiChaio/ccius f° r tne tidal radius is found to work well. Here 
r per i is the pericenter of the galaxy's orbit in the cluster, and cr\^\ 
and cr c i us are the velocity dispersions of the galaxy halo and the 
cluster, respectively. The tidal truncation of the dark matter halos 
of galaxies has been detected with combined s trong and weak 
lensing observations in several dense clusters ([N ataraian et all 
1998, 12002k iLimousin et ail 120071; lHalkola et all 120071) . Tidal 
radii of between 15 and 60 kpc have been found, in agreement 
with predictions. 

The case of M87 is not so simple though. M87 is at 
the center of at least a subcluster potential well, traced by 
the X-r ay emission and the dark matter mass profile de rived 
from it dNulsen & Bohringerl 1 19951: ISchindler et ail 1 19991 see 
Fig. [T2l . A galaxy at the center of its cluster experiences a 
symmetric gravitational field from the cluster d ark mat t er, an d 
is consequently not gravitationally truncated (Merritt, 1984). 
On the other hand, M87 has a relative motion of ~ 300 km/s 
with respect to the galaxies in the cluster core (Binggeli et al., 
1 19931) . and the galaxy distribution in the core is complicated 
and not centered on M87, containing a strong concentration 
around M84/M86 (Bing geli etail 119871) . It is possible that the 
M84/M86 concentration including the associated dark matter ex- 
erts a significant tidal field on M87. There is no obvious fea- 
ture in the density of the X-ray emitting gas at R tnmc ~ 150 
kpc around M87, but because the total mass within the trun- 
cation radius appears to be already dominated by cluster dark 
matter, a tidal truncation of the M87 mass distribution may be 
difficult to see in X-ray observations. However, if M87 was cur- 
rently tidally truncated by a tidal field with assumed mass center 
towards M84/M86, we would expect to see some of the tidally 
dissolved stars as PNe in our outer F7 fields at slightly redshifted 
velocities with respect to the systemic velocity. Within the lim- 
ited statistics, we do not see any PNs with such velocities out- 
side a projected radius of 161 kpc, see § 14.31 This suggests that 
if there has been a tidal truncation, it would have occurred some 
time ago during the interaction with another mass concentration. 
The most likely candidate in the Virgo core may be that around 
M84; at a relative velocity with respect to M87 of about 300 
km s" 1 , M84 could have travelled their current projected separa- 
tion in ~ 1 Gyr. 

On the other hand, due to the dynamical youth of the Virgo 
cluster, it is also possible that M87 has not been tidally af- 
fected yet, and is more similar to an isolated massive ellipti- 
cal galaxy. As already mentioned in § 14.31 the luminous halos 
aroun d isolated galaxies are expected from numerical simula- 
tions dAbadi et all 120061) to extend to the virial radius, i.e., to 
several hundred kpc, and well beyond their traditional luminous 
radius. Hence we now consider the possible origins of the trun- 
cation of the M87 stellar halo in the context of isolated galaxies. 

One possible explanation might lie in the fact that M87 is 
an old galaxy with a massive nuclear black hole, which points to 
much stronger AGN activity in the past than is apparent now. The 
feedback from its AGN through the surrounding hot gas might 
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at some redshift Zf have stopped the star formation in nearby 
satellite galaxies through, e.g., ram pressure stripping. When 
these satellites later accreted onto the galaxy, they would have 
predominantly added dark matter to the outer halo, so that the 
virial radius Ry of M87 kept growing, but the luminous radius 
stalled at Rv(z/)- On this assumption, we can estimate zj from 
the redshift dependence of the virial radius. 

The X- ray observations show that the hot gas extends ou t 
to 300 kpc dNulsen & BohringeA[l995tlMatsushita et alll2002l) . 
The derived integrated mass profile of the total gravitating mat- 
ter sh ows a change in slope at about 30 kpc (Matsushita et al.L 
120021) . and then increases linearly at large radii (see Fig. Q~2] 
above). The mass distribution inferred from the X-ray measure- 
ments thus provides evidence for two components: a galaxy dark 
matter compone nt and a cluster dark matter component. From 
the modeling of Nulse n~& Bohringerl (1 19951) and the rotation 
curve in Fig. [12] we can estimate the maximum circular veloc- 
ity generated by the M87 galaxy halo now to be v max ,M87 - 
400km s -1 . Using the results of lBullock et all d200ll) . this cor- 
responds to a present-day virial mass My(z = 0) 2.0 x 10 13 Mq 
and virial radius Ry(z = 0) - 470 kpc, several times larger than 
the truncation radius inferred from both the PN number counts 
and the outer halo dynamics, R lrunc - 149 kpc. For the same 
Vmax,M87, the virial radius of M87 would have been 149 kpc at 
redshift zj — 2.9, arguing that feedback would need to have been 
effective quite early-on to explain the observed truncation radius. 

A second possible explanation would assume that the accre- 
tion of dark matter and satellites onto M87 ceased with the col- 
lapse of the Virgo cluster core. In the new potential after the col- 
lapse, the satellites would both have been deflected from their 
nearly radial orbits with respect to M87, and have significantly 
larger impact velocities than previously, making accretion and 
merging with M87 suddenly less likely. The total mass of M87 
would thus not increase any further, stalling at the virial mass 
at that redshift. Moreover, the rotation curve in Fig. [12] shows 
that with the on-going collapse of the Virgo cluster a substan- 
tial cluster dark matter cusp has since built up within the halo of 
M87. The likely effect of this is an adiabatic contraction of the 
galaxy's luminous and dark halo. 

Th e two-component mass model of Nulsen & Bohr ingerl 
( 1995) (see equation[2]) for the present-day mass distribution pre- 
dicts within Rtrmc — 149 kpc, a galaxy mass of M(R tmnc ) M S7 - 
5.4 x 1O 12 M0 and a cluster dark matter mass of M(R trunc )yi rgo 
1 .2 X 10 13 Mq, assuming a flat rotation curve for the galac- 
tic contributio n to the mass. The lu minous mass of M87, (4- 
5)xl0 n M© dCappellari et all l2006h . is consistent with the es- 
timated total galaxy mass. As an example, consider truncation 
of the accretion onto M87 by the collapse of the Virgo core at 
redshift z = 0.5. Using a lower v c>mXj M87 a 300km s -1 for the 
M87 halo before ad iabatic contraction and again the results of 
iBullock et al.l d200ll) . the virial mass and radius at that redshift 
become M v (z = 0.5) 6.2 x 1O 12 M0 and virial radius R v (z 
0.5 ) =s 300 kpc. An (over)estimate of the adiabatic contraction 
can be obtained from angular momentum conservation at the 
outer radius, i.e., GM v (z = Q.5)R v (z = 0.5) = G[M{R tnm c)mi + 
M(R, runc ) vilso ]R, runc , giving R, runc /Rv(z = 0.5) = 0.3. This sug- 
gests that the observed truncation radius could well be the relic 
of the virial radius from the time when the cluster core collapsed. 

It is clear that more data are needed to pursue this question 
further. In particular, a larger number of PNs all around M87 
would be very useful to set stronger constraints on the tidal hy- 
pothesis. 



5. The luminosity-specific PN number for the M87 
halo and the ICL in Virgo 

The physical quantity which ties a PN population to the lumi- 
nosity of its parent stars, is the luminosity-specific PN number 
a = Npw/Lgai, where Npn is the number of all PNs in the popu- 
lation^] and Lg a i is the bolometric luminosity of the parent stellar 
population. Observations show that this quantity varies with th e 
(B - V) color of the stellar continuum light dHui et all 1 19931) . 
and simple stellar population models predict that it is a func- 
tion of the age and m etallicity of the parent stellar population 
(Buzzoni et all 120061) . Furthermore, within the framework of 
single stellar populations models, the a parame t er quantifies the 
avera ge PN lifetime t pn dVillaver et all 120021: ICiardullo etaTl 
2005) via the relation a = Btpn, where S is the 'specific evo- 
lutionary flux' and is nearly constant (see Buzz oni et all 2006, 
for a detailed discussion). The PN samples in elliptical galaxies 
and the ICL are all confined to the brightest 1 to 2.5 magnitudes 
of the PNLF Therefore we use the 0:2.5 parameter, defined in 
terms of N2.5, the number of PNs down to 2.5 magnitudes below 
the PNLF cut-off: 02.5 equals abo ut one tenth of a accor ding to 
the double exponential formula of lCiardullo et alj dl989l) for the 
PNLF. 

We can use the number of photometrically detected PNs in 
the FCJ, F7/F4 fields, and the luminosities of both the M87 halo 
and the ICL populations sampled in the surveyed areas, to com- 
pute the «2.5 values for these two components. Since the M87 
halo and the ICL coexist at the two field positions, we deter- 
mine the fraction of the PNs in the photometric sample bound to 
the M87 halo or in the ICL according to the fraction of spectro- 
scopically confirmed PNs associated with each component in the 
LOSVD, i.e., with the narrow M87 peak or the nearly uniform 
velocity distribution for the ICL. 

Luminosity of the ICL - We estimate the luminosities of the 
ICL stellar populat ion in our fields from the deep photometry of 
Miho s et ail d2005l) . Such surface brightness measurements for 
the diffuse light generally contain the cumulative contributions 
from extended galaxy halos, the true ICL, and from excess un- 
resolved background galaxies above that adopted from the sky 
subtraction. Depending on the method of sky subtraction, the 
homogeneous part of the last compon ent may be include d in the 
sky measurement. The photometry of Miho s etaLI d2005l) shows 
a 'plateau' of the surface brightness at a value of /uy = 27.7 
half way between M87 and M86, where the F7/F4 fields are sit- 
uated. As we have seen in Sections ^. 3114. 41 there is no contribu- 
tion from the halo of M87 to the plateau. IWilliams et all (2007) 
estimate the surface brightness of background galaxies from 
their deep imaging survey with the Hubble Space Telescope's 
Advanced Camera for Surveys (ACS), in a small (intracluster) 
field within our field F4. These galaxies, which are resolved in 
the ACS data, would contribute a diffuse surface brightness of 
fiy — 28.6 i n ground-based da ta. The sky subtraction procedure 
adopted by iMihos et all d2005l) would have already subtracted 
this component if its surface brightness is similar in the edges 
of the mosaic where the sky was measured. In the following, we 
therefore use fiy - 27.7 for the ICL surface brightness in this re- 
gion, with a possible uncertainty of a few tenths of a magnitude 
due to a possible inhomogeneity of the background sources. 

The ai.5 values and their implications - In what follows we 
make the assumption that the ICL surface brightness is constant 
in the FCJ, F7/F4 fields at this value of fiy ~ 27.7 mag arcsec" 2 . 



5 This is given by the integral over the whole eight magnitude range 
of the Planetary Nebula Luminosity Function (PNLF). 
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In Table [3] we give the corresponding ICL luminosities and the 
number of spectroscopically confirmed PNs from the ICL in 
these fields, N spectr . In Table [4] we give the number of PNs in 
the complete photometric samples in the overlap area covered 
by the spectroscopic follow-up, N p i w i, and the fraction of PNs in 
the M87 halo and the ICL according to their measured LOS VDs. 
Because of the small number statistics, we compute 0-2.5 for the 
M87 halo in FCJ, and for the ICL in F7, where their respective 
contributions are largest. 

Table 3. ICL PNs and sampled luminosities in the colored re- 
gions of Fig.fTTI 



Field 


N 

1 " spectr 


Area 


L 




kpc 2 


Ly,o 


FCJ red 


2 


4157 


1.2 x 10" 


F7 red + green 


7 


13258 


3.8 x 10" 



Uncertainties in the a-2.5 values - The luminosity of the M87 
halo is computed using Mont e Carlo integration of th e Sersic fit 
to the surface brightness from lKormendv et a l. (2008) in the FCJ 
field, and the errors here are of the order of few percent. 

The lum i nosity of the ICL is computed using fiv = 27.7 from 
Miho s et ail d2005i) . We independently estimated the ICL sur- 
face brightness by comparing the reprojected s urface brightness 
profile of the M87 halo with the Sersic fit of iKormendv et all 
(2008) in the F7 field. This results in an azimuthally averaged 
ICL jiy — 28.5. W hile 50% fainter than the measurement of 
Miho s et al] (l2005), the two values may be quite consistent when 
taking into account that the ICL is observed mostly on the side 
of M87 towards M86/M84. 

Considering the uncertainties in the surface brightness for 
the M87 halo and ICL, and the statistical errors in the number 
of detected PNs, the a 2 .5 values for M87 and ICL differ at the 
~ 2cr level. We speculate that a^.sjch is a factor 2 larger than 
#2.5^87 because of different metallicity distributions in the ICL 
and the M87 halo, with a larger fractio n of metal-poor stars in 
the intracluster component, as shown by Willia ms" et alj ([2007 ). 



Table 4. Number of PNs in the photometric samples, for the M87 
halo and the ICL. 



Field 












M87 


ICL 


FCJ 


16" 


14 


2 


F7 


24* 


10 


14 



" Number of PNs above the photometric completeness limit from 
Aguerri et al. (2005). Number of PNs above the photometric com- 
pleteness from TableQ]in F7_l and F7_2. 



The Nphot of each field is scaled by a factor 
$m" +2 ' 5 PNLF(m)dm 



PNLF(m)dm 



(6) 



where PNLF(m) is the analytic expression for the PNLF 
(Ciardull o et allfl989l) . M* and m* denote the absolute and ap- 
parent magnitude of its bright cutoff, respectively, and mu m is the 
photometric OT5007 limiting magnitude in each field. This scaling 
ensures that we account for all PNs within 2.5 mag of M* . 

The V-band luminosities in each field are converted to bolo- 
metric luminosities according to 



Lboi - £v©10 



-0.4(BC v +0.07) 



(7) 



According to lBuzzoni et al. (2006), a value of BCy = -0.85 mag 
can be taken as a representative correction for all galaxy types 
within 10% uncertainty. 

Finally, we obtain the bolometric luminosity-specific PN 
number 0-2.5: For the M87 halo light at the FCJ position it is 
«2.5,M87 = 3.1 x 10~ 9 PN L \ and for the ICL at the F7 position 
it is a-i.sjcL — 7.2 x 10~ 9 PN L^ 1 . The values of 0:2.5 for dif- 
feren t stellar populat i ons ar e well documen t ed dCiardullo et all 
120051; iBuzzoni et all 120061: ICoccato et all |2009j): a 2 . 5 in the 
range 3 - 10 x 10~ 9 PN L ' are observed for bright ellipticals 
and SOs. Both the 0-2.5 values obtained for the ICL and for the 
M87 halo stars are consistent with those of old (> 10 Gyr) stel- 
lar populations. 



6. ICL and the dynamics of the Virgo cluster core 

6. 1 . ICPNs and dwarf spheroidals 

The LOSVDs in Figure [6] show the dynamical components in 
the Virgo cluster core: the halo of M87, and the ICL component 
traced by a broad PN velocity distribution. This component cov- 
ers the velocity range from 1 300 km s~ 1 down to the systemic ve- 
locity of M86 at -244 km s . Overlaying the spatial coordinates 
of the PNs on the d eep image of the Virgo cluster core (Figure|7l 
Mih os et ail d2005)) we can easily see the association of the PN 
components identified in the velocity - position space, with the 
morphological components of the surface brightness distribution 
in the Virgo core. The M87 PNs are confined to still relatively 
bright regions covered by the M87 halo, while the ICL PNs are 
scattered across the whole region. 

For comparison, we examine the phase space distribution of 
dwarf elliptical galaxies in a region covering our target fields 
for the PNs spectra, i.e., in a 1.5 degree diameter circle cen- 
tred on the midpoint of M87 and M86 (Figure |7j. The aim is 
to search for possible associations between our ICPNs and the 
positions and velocities of the dwarf galaxies. The top panel of 
Figure [6] shows a histogram of the LOS velocities for all dwarfs 
in the region marked in Figure|7] The velocities form a flat, uni- 
form distribution extending to larger positive velocities than the 
ICPNs. We ask whether any of the PNs could be physically as- 
sociated with the dwarf galaxies. There are only two potential 
associations: one is between a PN with velocity 818km s -1 and 
a close-by dE at 791km s _1 . This dE has total blue apparent 
magnitude 15.4, i.e., M#=-15.4 8 (using the assumed di stance) or 
L B =IA x 10 8 L Q . According to lCiardullo et al.1 f2005). galaxies 
fainter than M# ~ -19 and bluer than V - 1 < 1.1 produce about 
one [O in] -bright PN in every 4 x 10 8 Lq ■ It is therefore unlikely 
that the dwarf galaxy in question produced the PN detected here 
(as the expected number is 0.35) although it cannot be ruled out 
completely. On the other hand, the second association may well 
be genuine: this is of a PN with velocity +28km s -1 and a close- 
by Sb spiral galaxy to its west, which has velocity 30km s _1 and 
fi, o ,= 10.91 (and therefore is capable of producing between 85 
and 2000 PNs depending on its age, see Buzzoni et al. 2006). 

Although the majority of the ICPNs do not appear to be 
physically associated with the dwarf galaxies (i.e., the PNs are 
unlikely to originate in the dwarfs), their distribution in veloc- 
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ity space is at least partially similar, indicating that they follow 
similar dynamics. 

6.2. Dynamical status of the Virgo core 

The velocity distribution of dwarf spheroidals (dE+dSO) in a 2° 
radius circular region centred on M87 is very flat and broad, with 
the peak of the dist ribution at 1300 km s" 1 and a long tail of 
negative velocities (Bi nggeli etafl Il993l) . The LOSVD of the 
ICPNs now confirms that this asymmetry is also present in the 
very center of the Virgo core, in a region of 1° diameter. Fig. [8] 
shows that velocities near the systemic velocity of M86 are seen 
to about half-way from M86 to M87. 

The asymmetry and skewness of the LOSVD may arise 
from the merging of subcluste rs along the LOS as described by 
ISchindler~ & Bohringer ( 1993). In their simulations of two merg- 
ing clusters of unequal mass, the LOSVD is found to be highly 
asymmetric with a long tail on one side and a cut-off on the other 
side, shortly (~ 10 9 yrs) before the subcluster s merge. 

The observed LOSVDs of the PNs, GCs (ICote et all l200lh . 
and (dE+dSO) in the Virgo core may therefore be interpreted as 
additional evidence that the two massive subclusters in the Virgo 
core associated with the giant ellipticals M87 and M86 are cur- 
rently falling towards each other - more or less along the LOS, 
with M87 falling backwards from the front and M86 forwards 
from the back - and will eventually merge, i.e. the entire core 
of the Virgo cluster must then be out of virial equilibrium and 
dynamically evolving. 

The distribution of the brightest galaxies in Virgo also favors 
a recent and on-going assembly: West & Blakeslee (2000) found 
that Virgo's brightest elliptical galaxies tend to be aligned along 
the principal axis of the cluster (which is inclined by only about 
10-15 degrees to the line of sight) and which on larger scales 
connects Virgo to the rich cluster Abell 1367. This work suggests 
that the formation of the cluster is driven by infall along this 
filament. 

Do the halos of M87 and M86 already touch each other, 
or are they just b efore their close pass? PNLF distances 
dJacobv et al.Lll990|) and ground-ba sed surface brightness fluctu- 
ation distances dTonrvet all 120011) indicate that M86 is behind 
M87 by just under ~0.15 magnitudes. The globular cluster LF 
turnover also suggests that M86 is likely 0. 1 to 0.2 magnitudes 
more distant than the main body of Virgo dKundu & Whitmora, 
120011) . However, the most recent surface brightness fluctuation 
measurements by Mei et al. (2007) find that M87 and M86 are 
only at very slightly different distances. Within the errors, the 
distance moduli (M87: 31.18 + 0.07, M86: 31.13 + 0.07) are 
consistent with being either at the same distance or separated by 
1-2 Mpc. Unfortunately the evidence from the relative distances 
of M87/M86 is not conclusive at this stage. 

6.3. Implications for the formation of the ICL 

The observational facts concerning the ICL in the Virgo cluster 
core are: 

1. The LOSVD of the ICPNs is not symetrically distributed 
around the systemic velocity of M87. Those between M87 
and M86 are mostly at 'bluer' velocities, i.e., < 800 km s . 
"Red" velocities are only seen in the field 0.°8 north of the 
line connecting M87 with M86; see Fig. [7] 

2. While the dwarf spheroidals' (dE+dSO) LOSVD in the re- 
gion marked in Figure [7] extend into 'red' velocities, up 
to 2500 km s" 1 , ICPNs with velocities greater than 1800 



km s are seen only at its northern perimeter, while those 
in the region between M87 and M86 are confined to < 800 
km s~'(see Figs. [6] |7] EJ. This is not a consequence of 
the filter used in the photometric selection of these objects, 
which still ha s a transmission of 50% to [OIII]/15007 at 
~2275kms- 1 dFeldmeier et al.U2003l) .Fl 

3. The morphology of the ICL between M87 and M86 is 'dif- 
fuse'; it is mostly not in tidal tails or streams ( Miho set all 
120051) . 

4. The measured C2.5 parameter for the ICL is in the observed 
range for old stellar populations. 

5. The metallicity distribution of t he RG B stars associated with 
the ICL in the IWilliams etail d2007l) field i s broad, with a 
peak a t about 0. 1 solar, and the best model of lWilliams et al.l 
(2007) indicates an old stellar population (> lOGyr). 

Point 1) indicates that the ICL did not have the time to phase 
mix yet, supporting the idea that much of this diffuse component 
is falling towards M87 together with the M86 group. Points 2) 
and 3) argue against the origin of this ICL from cu rrent harrass- 
ment of dwarf irregulars by the cluster potential dMoore et all 
11996b . firstly because none of the ICPNs between M87 and M86 
are found at velocities redder than 800 km s _1 , while we see 
nearby dwarf galaxies in the 1800 < Mlos < 2500 km s _1 range, 
and secondly because the diffuse morphology of the ICL sug- 
gests that it has not recently been dissolved from the dwarfs, but 
had time to phase mix at earlier times in the sub-component of 
the Virgo cluster which is now falling towards M87. Points 4) 
and 5) indicate that the parent stellar population of the ICPNs 
is an old population, and point 5) makes dwarf sp heroidals un- 
likely progenitors (recall that the dSph analyzed bv lDurrell et al.l 
(120071) had a narrow metallicity distribution at [Fe/H]^ -2). 
Dwarf elliptical galaxies have a wider metallicity distribution, 
and they could be disrupted during passages through the high- 
density regions around M87. Some of the observed ICPNs could 
trace stars on a stream from such a disruption event (particularly, 
those at Vlos ~ 800 km s _1 ), but those at low velocities in Fig. [8] 
are falling in with M86, so have not crossed the dense regions 
yet. Point 3) then implies that these must have been part of the 
diffuse halo of the M86 group. 

We conclude that we have found observationa l evidence in 
the Vi rgo core for the mechanism described by Rudic ket al.l 
(2006): we observe the diffuse component 'pre-processed' in 
the M86 sub-group, which is or has been gravitationally un- 
bound from M86 as this substructure is being accreted by 
M87. The idea that the diffuse light is being stripped from 
the M86 sub-group is consistent with the observed highly 
sk ewed LOSVD and with the pr edictions from the simulations 
of lSchindler & Bohringer! (119931) . Note that the light in the M86 
subgroup is tidally stripped by the more massive M87 compo- 
nent, while these two sub-structures merge along the LOS; we 
do not see a diffuse ICL with a broad velocity component red- 
wards of the systemic velocity of M87, because it has not yet 
been formed. 

This scenario is also consistent with the simulations of 
Mur ante et al.ld2007l) . Their statistical analysis of the diffuse star 
particles in a hydrodynamical cosmological simulation indicates 
that most of the ICL is associated with the merging tree of the 
brightest cluster galaxy, and about 80% of the ICL is liberated 
shortly before, during and shortly after major mergers of mas- 

6 There may be a different, small selection effect due to the finite 
limiting magnitude of the photometric PN survey, as the ICPN that we 
detect are slightly biased towards objects on the near side of the cluster 
core. 
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sive galaxies. The results from M urante et al.l d2007l) imply that 
the main contribution to the ICL comes from merging in earlier 
sub-units whose merger remnants late r merge with the final cD 
galaxy. Similarly, Rudic keit al.l d200 6) predict that violent merg- 
ing events quickly add ICL, and without or between these events, 
the ICL fraction rises only slowly. Once the M86 subgroup has 
finally merged with M87, this will have created the most mas- 
sive galaxy in the then Virgo cluster, and the ICL in the future 
Virgo core will indeed have originated mainly from the progen- 
itors associated with its merging tree. 

7. Summary and Conclusions 

Using high resolution multi-object spectroscopy with 
FLAMES/MEDUSA on the VLT we confirm a further 12 
PNs in the Virgo cluster, located between 130 and 250 kpc from 
the center of M87, and obtain their radial velocities. For most 
of these objects we also detect the second line [Om]/l 4959A. 
These PNs trace the kinematics of diffuse light in Virgo, at 
typical surface brightness of yUy = 27.5. 

The phase-space distribution for the new sample of PN ve- 
locities combined with earlier measurements at 60 and 150 kpc 
from M87 illustrates the hierarchical nature of structure forma- 
tion. One group of PNs has an unrelaxed distribution of veloc- 
ities with a range characteristic for the still assembling Virgo 
cluster core, while the second group has a narrow velocity dis- 
tribution which traces the bound, cold outer halo of M87. We 
summarize our results for these two groups in turn. 

7.1. Dynamical status of the Virgo cluster and origin of the 
ICL 

Seven of the newly confirmed PNs are genuine intracluster PNs 
in the Virgo core, not bound to M87. Their spatial and veloc- 
ity distribution indicates that we are witnessing the gravitational 
stripping of the diffuse light component around the M86 group, 
as this sub-structure is being accreted by the more massive M87. 
We do not see a diffuse ICL with a broad velocity distribution 
including red-shifted velocities around M87, because it presum- 
ably has not been formed yet. 

On the basis of the LOSVDs of ICPNs and galaxies in the 
Virgo core, we surmise that M87 and M86 are falling towards 
each other nearly along the line of sight, and that we may be 
observing them in the phase just before the first close pass. We 
thus conclude that the heart of the Virgo cluster is still far from 
equilibrium. 

Finally, the QT2.5 values determined for the ICL indicate an 
old stellar population . This is consistent with the analysis by 
IWilliams et al.l d2007l) of the colour-magnitude diagram (CMD) 
for ICL red giant stars, which showed an old stellar population 
(> 10 Gyr) with a large spread of metallicities. Differently, the 
CMD for a nearby dwarf spheroidal gala xy indicates a simi - 
larly old, but metal-poor stellar population dDurrell et al.L 12007). 
Together with the observed Mlos distribution of the ICPNs, these 
results suggest that at least some of the ICL in Virgo originates 
from stars unbound from the brightest and most massive galax- 
ies. 

7.2. The M87 halo 

The other five of the newly confirmed PNs are associated with 
the bound halo of M87, at a mean projected radius R = 144 kpc 
from the centre of the galaxy. These PNs have velocities close to 



the systemic velocity of M87, with a small dispersion, and are 
furthermore segregated spatially from the rest of the (intraclus- 
ter) PNs, as shown in Figs.|7]and[8] 

The LOSVDs of the M87 PNs both in the new fields and in 
the FCJ field of A04 are consistent with no rotation of the outer 
halo around the photometric minor axis, outside R = 15 kpc. We 
cannot test whether the halo is rotating around the major axis. 
The rms velocity dispersion for the 5 M87 PNs at 144 kpc is 
78+25 km s , much smaller than the central velocity dispersion. 
Together with the results of A04, this indicates that the M87 halo 
becomes 'dynamically' cold beyond 50 kpc radius. 

The PNs around the systemic velocity of M87 are confined 
to radii R < 161 kpc. The absence of M87 PNs at larger radii 
with respect to th e extrapolated Sersic fit to the surface bright- 
ness profile from Kor mendv et alj d2008l) . despite being based 
on small numbers, is significant at the ~ 2cr level. This suggests 
that the edge of M87 has been detected, and it occurs at quite a 
large average ellipse radius - R lnmc = 149 kpc. 

We have tested the hypothesis of a truncated stellar halo dy- 
namically, using the observed stellar kinematics of the M87 halo. 
Combining all the velocity dispersion data available in the liter- 
ature with our new M87 cr PN data, we have solved the spherical 
Jeans equation assuming a total gravitational potential as traced 
by the X-ray emission. Within this framework, the Jeans model 
is able to reproduce the 'cold' PN velocity dispersion in the outer 
halo of M87 only if (i) the orbital structure in the outer halo be- 
comes highly radially anisotropic, with f3 > 0.4 at r > 10 kpc, 
and (ii) the intrinsic luminous density is truncated. This dynami- 
cal argument can be circumvented only if the stellar distribution 
were strongly flattened along the line-of-sight in the surveyed 
fields. At these radii this flattening would be local and imply 
that the spheroidal stellar halo would end at even smaller radii. 

The evidence for the truncation of the luminous halo of M87 
thus comes from both the spatial distribution of the PNs with 
velocities near the systemic velocity of M87 (Fig. |8j, and from 
the small velocity dispersions in both the A04 field and in the 
new fields near the outer edge. 

The reason for the truncation is not obvious; we discuss pos- 
sible mechanisms in Section [431 Differently from some dense 
clusters where lensing analysis indicates that galaxies outside 
the cluster center are tidally truncated by the dark matter cusp of 
the cluster, M87 in Virgo is located at the center of the deepest 
potential well traced by the X-ray isophotes. We also do not see 
unbound PNs with velocities near 1300 km/s further out from 
M87. This suggests that if there has been a tidal truncation, it 
would have occurred some time ago during the interaction with 
another mass concentration, such as around the other massive 
galaxy in the Virgo core, M84. Alternatively, due to the dynam- 
ical youth of the Virgo cluster, it is also possible that M87 has 
not been tidally disturbed yet, and is more similar to an isolated 
mass ive elliptical gala xy and should thus still be accreting mat- 
ter dAbadi et al.Ll2006l) . In this case possible explanations for the 
truncation could be early AGN feedback effects that indirectly 
truncate star formation in accreting satellites, or adiabatic con- 
traction of the M87 halo due to cluster dark matter collapsing 
onto the galaxy. 

The existing data cannot discriminate between these scenar- 
ios. The next step in this project is therefore to obtain a sample 
of at least a few hundred measured PN velocities covering the 
whole M87 halo. This is required to verify that the velocity dis- 
persion decreases everywhere around the galaxy, and to obtain 
statistically better constraints on phase-space structures in the 
surrounding ICL, including possible stars tidally dissolved from 
M87. With a homogenous imaging and spectroscopic PN survey 
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covering the whole halo of M87 out to 40 arcmin we will be 
able to accurately measure the rotation, radial anisotropy of the 
orbits, and truncation of the outer halo of M87. 
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